Structure of divergences in Drell-Yan process 



with small transverse momentum 

Junegone Chay 1 ']^ and Chul Kim 2, [T] 

1 Department of Physics, Korea University, Seoul 136-713, Korea 
2 Institute of Convergence Fundamental Studies & School of Liberal Arts, 
Seoul National University of Science and Technology, Seoul 139-743, Korea 

Abstract 

We consider the structure of divergences in Drell-Yan process with small transverse momentum. 
The factorization proof is not trivial because various kinds of divergences are intertwined in the 
collinear and soft parts at high orders. We prescribe a method to disentangle the divergences 
in the framework of the soft-collinear effective theory. The rapidity divergence is handled by 
introducing the 5 regulator in the collinear Wilson lines. The collinear part, which consists of the 
transverse-momentum-dependent parton distribution function (TMDPDF), is free of the rapidity 
divergence after the soft zero-bin subtraction. There still remains the problem of mixing between 
the ultraviolet and infrared divergences, which forbids the renormalization group description. We 
show that the mixing is cancelled by the soft function. This suggests that the collinear and 
soft parts should be treated as a whole in constructing a consistent factorization theorem. The 
renormalization group behavior of the combined collinear and soft parts is presented explicitly at 
one loop. We also show that the integrated PDF can be obtained by integrating the TMDPDF 
over the transverse momentum. 
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I. INTRODUCTION 



Theoretical predictions of high-energy scattering rely on the factorization of scattering 
cross sections, in which the hard, collinear and soft parts are separated to all orders in 
perturbation theory. Many inclusive cross sections have been proved to be factorized, and 
are used to compare with experiments. Less inclusive scattering processes are also of in- 
terest. Recently the transverse momentum distribution in the Higgs, Z boson, and it pair 
production has redrawn new interest both theoretically [IHS] and experimentally [61 [7] . But 
in high-energy scattering processes with nonzero momentum transverse to the beam direc- 
tion, the factorization proof becomes more involved. For example, in the Higgs production 
with the transverse momentum qx with Aqcd qr M, where M is the Higgs mass, the 
scattering cross section involves unintegrated, or transverse-momentum-dependent parton 
distribution functions (TMDPDFs). The study of the transverse momentum distribution 
was first heralded by Ref. [8]. However, the factorization proof, which appears ubiquitously 
in high-energy processes including Drell-Yan process, has been recently under debate again. 

Before we discuss many issues in proving factorization of the processes with transverse 
momentum, it is worth delineating why the Drell-Yan process with small transverse mo- 
mentum is so intricate. When we consider radiative corrections in massless gauge theories, 
there appear ultraviolet (UV) and infrared (IR) divergences. The UV and IR divergences 
should appear separately with no mixing in order for the theory to be consistent. In many 
inclusive processes, this was shown explicitly in the collinear and soft parts. In the process 
where the transverse momentum for a final-state particle is fixed, the UV divergence to be 
obtained by integrating over the transverse momentum does not emerge yet. Therefore the 
separation of the UV and IR divergences seems murky at first sight. Furthermore, when 
the transverse momentum is fixed, there appears another type of divergence called rapidity 
divergence, or lightcone singularity. The issues in the Drell-Yan processes with small trans- 
verse momentum reside in the appropriate treatment of rapidity divergence and in finding 
an adequate formalism in which the separation of the UV and IR divergences is established. 

The first issue is the existence of rapidity divergence, which has been known for a long time 
[9] , and various ways of handling it have been suggested. The source of the rapidity diver- 
gence can be understood by the following argument: In QCD, the Kinoshita-Lee-Nauenberg 
theorem (TU1 E] guarantees that the divergence in virtual corrections at a given order in 
the strong coupling a s is cancelled by the divergence in real gluon emissions, thus rendering 
inclusive scattering cross sections free of infrared divergence. The presence of infrared diver- 
gences can be traced by considering the kinematics in each case. In virtual corrections, the 
loop momentum can be collinear to an energetic particle, or soft, which results in collinear 
or soft divergences, or it can be both collinear and soft. In real emissions, the emitted gluon 
can also be collinear to an energetic particle or soft, or both as long as the phase space allows 
it. These divergences cancel in inclusive processes, and the inclusive scattering cross sections 
are physically meaningful in the absence of the collinear and soft divergences. However, if 
we consider differential cross sections with fixed transverse momentum, this cancellation 
becomes incomplete. That is, when the transverse momentum of an emitted gluon is fixed, 
the phase space of the emitted gluon does not cover all the available phase space in virtual 
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corrections, and the IR divergences do not cancel because of the imbalance of the phase 
space between real and virtual corrections. The presence of the infrared divergence due to 
this incomplete cancellation is referred to as the "rapidity divergence" or lightcone singu- 
larity since it is the result caused by the interaction of the energetic particle with collinear 
gluons with infinite rapidity. The dimensional regularization cannot regulate the rapidity 
divergence. 

The difficulty in regulating the rapidity divergence was first noticed by Collins [9 J . Collins 
and Soper [12l [13] suggested an intriguing idea of defining the operator whose matrix el- 
ements generate TMDPDF by tilting the Wilson lines which connect the fermion fields to 
make the operator gauge invariant. In association with this tilting, they introduced an addi- 
tional scale of which the TMDPDF should be independent, and the corresponding evolution 
equation was developed. This method, while effective in handling the rapidity divergence, 
has a nuisance. First, the dependence of the additional scale appears in the hard, collinear 
and soft parts, which troubles the factorization proof. Though the dependence of the addi- 
tional scale disappears when all the parts are summed, it is not clear how to approach the 
physical limit, that is, the limit in which all the particles are put on the lightcone. Finally 
it is not straightforward to obtain the integrated PDF from the TMDPDF by integrating 
over the transverse momentum [T3] . 

Chiu et al. [151 IB] developed the idea of rapidity renormalization to handle the rapidity 
divergence. In high-energy scattering, there are collinear particles and soft particles, and 
they note that particles with the same offshellness can be labeled as collinear or soft depend- 
ing on their rapidities. Therefore a rapidity scale is set up to distinguish collinear and soft 
particles. However, since the physics is independent of this arbitrary scale for the separation, 
the evolution equations of the collinear and the soft parts can be obtained with respect to 
the scaling of this rapidity scale. Technically they modified collinear and soft Wilson lines 
in such a way that Wilson lines depend on the rapidity scale. The rapidity scale dependence 
of the Wilson line is determined to satisfy the requirement that the cross section be inde- 
pendent of the rapidity scale, that is, the rapidity scale dependence of the collinear part is 
cancelled by that of the soft part. Though the motivations come from different reasoning, 
the approaches of Collins et al. [5] and Chiu et al. [T5J [TH] share in common the fact that 
they introduced an additional scale of which physical results should be independent. This 
is reflected in the evolution equations of the collinear and soft parts. The purpose of this 
approach is also focused to handle the rapidity divergence. 

In this paper, a straightforward way of regulating the rapidity divergence is presented in 
the framework of the soft-collinear effective theory (SCET) [T7HT9] . The regularization for 
the rapidity divergence is achieved by introducing the 5 regulator in the collinear Wilson 
lines only and not in any other propagators. Chiu et al. [20] introduced the 5 regulator, 
which is an infrared regulator inserted in every propagator, as well as in the Wilson lines. 
They showed that the 5 dependence of the collinear part in the back-to-back current is 
cancelled by the zero-bin subtraction. They stressed the importance of the proper zero- 
bin subtraction in the cancellation of the 5 dependence. The 5 regulators can be given as 
arbitrary parameters, or they can be related to the offshellness of the external particles. 
In either case, the 5 dependence in the collinear and soft parts is different and cannot be 
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cancelled unless a specific relation is imposed. A simple relation can be found in a special 
case like a back-to-back current when the offshellness is employed to regulate divergences. 
However, when multijets are involved, 5 depends on the offshellness of all the external 
particles in a complicated way. The collinear Wilson line in the n direction is obtained by 
considering the emission of n-collinear gluons from the other particles not in the n direction. 
Then the intermediate states become offshell, and the collinear Wilson line is obtained by 
integrating out this offshellness, and taking the leading term. If a back-to-back current 
is involved only, the 5 regulator can be related to the offshellness of the particle in the 
n direction for the n-collinear Wilson line. If we consider a process in which there are 
many jets, the offshellness in a collinear Wilson line depends on all the other jets in the 
process. On the other hand, the n-soft Wilson line is obtained by considering soft gluon 
emissions from the n-collinear particle. Therefore if the offshellness is included, the soft 
Wilson line involves only the offshellness of the corresponding collinear particle. Here the 
rapidity divergence is handled by inserting the 5 regulator in the collinear Wilson lines only. 
Other regulators, though similar to the 5 regulator in the collinear Wilson lines in form, 
regulate IR divergences. 

The second issue is to separate the remaining UV and IR divergences. The presence of the 
mixing between the UV and IR divergences is troublesome because the physics at a high scale 
and a low scale is mixed, which does not make sense. Therefore the decoupling of the UV and 
IR divergences in loop calculations is essential to guarantee the consistency of the theory. In 
contrast to inclusive processes, the appropriate combination in which this decoupling occurs 
is still vague. It makes the definition of the TMDPDF itself controversial. The problem is 
whether the soft contribution or part of it should be included in the matrix element of the 
collinear operators which comprise the backbone of the TMDPDF. The underlying idea is 
to find a suitable combination of collinear and soft parts such that UV and IR divergences 
at higher orders are separated in order to prove factorization. A specific combination of 
collinear and soft parts has been claimed to achieve this separation [2Tj . But once soft part 
is included in the TMDPDF, there follows another issue, the universality of the TMDPDF. 
The soft part interacts with all the different collinear sectors, and it is different for different 
final hadronic states. Therefore the definition of the TMDPDF in one process may lose 
its meaning in other processes once the soft part in a specific process is included in the 
TMDPDF though the TMDPDFs in two different processes can be related to each other. 
This kind of relation for the PDF in deep inelastic scattering and Drell-Yan processes has 
been discussed when the soft part is included in the PDF [22]. 

In Ref. [23], an intriguing method to define the factorized TMDPDF was developed. First, 
the authors put all the particles on the lightcone, which is physical. Therefore there is no 
need to introduce additional scales due to the tilting of Wilson lines or rapidity separation. 
Obviously a technique to handle the rapidity divergence is necessary, and the 6 regulator, 
slightly different from that in Ref. |20j , is employed. They defined the TMDPDF in terms of 
the collinear part with the soft function, which is similar to the approach of Collins |21| . The 
reason is that the 8 dependence, which indicates the presence of the rapidity divergence, is 
cancelled only through a specific combination of the collinear and soft parts as they claimed. 
It is based on the observation [24H26] that the zero-bin contribution can be identified as the 



4 



soft part. It is true for back-to-back currents to set up a relation of the regulators in the 
collinear and soft parts, but in general it does not hold. For example, if there are many jets 
or heavy colored particles in the process, there is no simple relation between them. 

In this paper, the TMDPDF is defined in terms of the collinear fields only, thus it is 
universal in all high-energy processes. All the particles are put on the lightcone and the 
rapidity divergence is handled by the 8 regulator in the collinear Wilson lines. The collinear 
part itself along with the corresponding zero-bin subtraction is independent of the 8 regu- 
lator. The zero-bin contribution is performed in the limit where collinear particles become 
soft (not ultrasoft) to avoid double counting [27] . The infrared divergence from the soft 
Wilson lines is cancelled by that of the collinear part, which is controlled by the offshellness 
of the external particles. Combining all the ingredients, we obtain a factorized form of the 
scattering cross section with small transverse momentum in Drell-Yan processes. The hard 
part is the Wilson coefficients obtained in matching the current operators between QCD 
and SCET. The collinear part consists of the TMDPDF, and the remaining part is the soft 
function. We develop a method to handle the rapidity divergence, with the decoupling of 
the UV and IR divergences, and show the result explicitly at next-to-leading order (NLO) 
accuracy. 

The structure of the paper is as follows: In Section [TTJ the factorized form of the scattering 
cross section in Drell-Yan process is presented. In Section III, the TMDPDF is defined as 
the matrix elements of collinear fields, and the one-loop corrections are computed. The 
role of the zero-bin subtraction to avoid double counting is discussed and we argue that 
the soft zero-bin subtraction is the appropriate method for small transverse momentum 
qr ~ QX. In Section IV, the soft function is defined in terms of the soft Wilson lines, and 
its one-loop correction is presented. The relation between the TMDPDF and the integrated 
PDF is discussed in Section [V] and we explain how the integrated PDF is obtained from 
the TMDPDF. In Section |VI[ the renormalization group behavior of the scattering cross 



section at NLO is presented. In Section [VTT] the decoupling of the UV and IR divergences is 
elucidated in detail. And finally in Section VIII, we summarize the procedure for taming the 
rapidity divergence with the decoupling of the UV and IR divergences, and give a conclusion. 
In Appendix, the properties of the /i 2 distribution and the infinity distribution are described. 



II. FACTORIZATION IN DRELL-YAN PROCESS 

Let us consider Drell-Yan process pp — > £ + £~ + X, where the lepton pair £ + £~ is produced 
with small transverse momentum q T (Aq CD <C q T <C Q). The invariant mass of the lepton 
pair Q is the large scale. The incoming particles are n-collinear and n-collinear, and their 
momenta scale as 

Pn = (n- Pn,n- p n ,Pn±) ~ Q(l, A 2 , A), p% = (u • p n , Tl • p n , PnA-) ~ Q(X 2 , 1, A), (1) 

where A = qr/Q is the small parameter. We also denote n ■ p = p~ , n ■ p = p + . There 
are soft particles, whose momenta scale as p^ ~ Q(X, A, A), in the final state to balance the 
transverse momentum of the lepton pair. In SCET, the interactions between collinear and 
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soft particles are not allowed since they put each particle far off their mass shells. Therefore 
the soft part is decoupled and is expressed in terms of the soft Wilson lines. The ultrasoft 
(usoft) particles can also exist, but they do not contribute to the cross section with transverse 
momentum q T ~ QX. 

The differential scattering cross section can be written as 



da 



s (2%y 



d^xe-^i-g^N^NzljUxpviO^Nz), 



(2) 



where o"o = 4ira 2 Q 2 /(3q 2 N c ), with the electric charge Qf and the number of colors N c . The 
electromagnetic current is given by J /t = fj^f. The sum over the flavors / is implied. We 
consider only the electroproduction here, but the weak interaction can be easily implemented. 
We choose the direction of the momentum Pi for the nucleon Ni to be in the n direction, and 
that of P 2 for N 2 in the n direction. In SCET, the current operator for the quark- ant iquark 
annihilation can be written as 



J„(x) = C(Q)e- J ^^-/ 2+ "^'^ 2 )x„,„>; t 7M^Xn,p(x ± ), 



(3) 



where Q 2 = q 2 , and Xn = is the gauge- invariant collinear fermion field with the 

collinear Wilson line W n . C(Q) is the Wilson coefficient obtained in matching between the 
full theory and SCET. Unlike the conventional SCET formalism in which the momenta of 
order Q and QX are label momenta, only the momenta of order Q are the only label momenta 
to be extracted in Eq. ([3]). It is because the soft momentum of order QX is the dynamical 
degree of freedom, and the size of the fluctuation after extracting label momenta is of order 
x± ~ (QX)- 1 . 

The differential cross section in SCET can be written as 



da 



d A q o"o 

(2~ 



7TI 4 S 



-gTN c ) / dutdu 2 / d 4 xe~ tq - x e 



-iq-x i(uJin-x/2+uj2n-x/2) 



\C(Q)\ : 



x (N 1 (P 1 )N 2 {P 2 )\[xJ(cj 1 - n ■ V r )]Y^ ± ,Y n [6(u 2 + n ■ V) X n(*±) 
xx n Yh±uY nXn (0)\Ni(Pi)N 2 (P 2 )), 



(4) 



where lo\ = x^n ■ Pi, u 2 = x 2 n ■ P 2 , and Xi, x 2 are the longitudinal momentum fractions of 
the incoming partons. Here n-V, n-V are the label momentum operators for the particles in 
the n and n directions respectively and the operators are applied inside the brackets. Since 
the soft Wilson lines are decoupled from collinear particles, they can be extracted out, and 
expressed in terms of the vacuum expectation values. The fields at xi can be expressed by 
the fields at the origin as 



d 2 k 1± e~ ikl ^ x 
d 2 k 2 ±e 



2i, „-ik 2 ±-xx 



S {2 \k 2± +V A )(xn)y) 



(Y^Y w ) ab (^ ± )(Y n Y n ) cd (0) = J d 2 r 1 ^^(Y} L Y n ) ab 5^ 2 \r 1± +iV i _){Y^Y n ) cd ^), (5) 
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where 'Pi is the operator extracting the transverse momentum. The TMDPDFs are defined 
after taking the spin average as 

(Nl\(Xn) a a [(XnP( X l ~ j^)^^ ~ W = ^aT^ (f) ^ ' ^ ' 



(N 2 \[s(x 2 + ^)5 (2) (k 2 ± + Pjfe)«]fe)J|iV 2 ) = r ^^ 6 (|) a /f/iv 2 (^2,k 2± ). (6) 



After straightforward algebra, the differential scattering cross section can be written as 
da 

x fq/N 1 (xi,ki±)fq/N 2 (x 2 ,'k2±)S(rij), (7) 



J dx 1 dx 2 cr H(Q 2 ,n) J d 2 k 1± d 2 k 2± d 2 ri ± 6 ( - 2 \k lx + k 2± - rj ± - q ± ) 



where s = n ■ P\n ■ P 2 , Q 2 = XiX 2 s, and qj_ is the transverse momentum of the lepton pair. 
From Eq. @, the TMDPDFs are given as 

f q/Nl (x, k ± ) = (N 1 (P 1 )\x n ^ 2 \k ± - V^5{n -V-xn- P 1 ) Xn \N 1 (P 1 )), 
f- q/N2 (x, kj.) = (A^ 2 (P 2 )|xJ^ (2) (k± + 7>j)<y(n • V + x« • P 2 )x^|iV 2 (P 2 )), (8) 
and the soft Wilson line is given by 

SM = ^(o|tr[r n ty^ 2 W+ iV±)YiY n ] |0). (9) 



c 



There is an alternative definition of TMDPDF in previous literature in terms of the delta 
functions of the transverse momentum in D — 2 dimensions instead of two dimensions. With 
this definition, the radiative corrections should change accordingly, but equivalent results can 
be obtained. However, the definition in Eq. ^ is preferred in this paper since the TMDPDF 
is an observable and it should be defined in four spacetime dimensions. If the Wilson lines Y n 
and Y n are usoft, the momentum fluctuation is of order QX 2 and the momentum operator in 
the delta function in Eq. ^ can be put to zero at leading order, and the usoft Wilson lines 
cancel. The fact that usoft interactions do not contribute to the process with the transverse 
momentum of order QX to all orders in a s manifests itself due to the cancellation of the 
usoft Wilson lines. 

Eq. ([7]) is the factorized form for the differential scattering cross section. The hard 
function H{Q) = \C(Q)\ 2 is the Wilson coefficient of the current, the collinear part consists 
of the product of two TMDPDFs, and the soft Wilson lines comprise the soft part. The 
radiative corrections of the hard, collinear and soft parts can be computed separately in 
perturbation theory. However, the infrared divergences, or the dependence of the regulators 
in actual calculations are intertwined between the collinear and soft parts, of which the 
disentanglement is the main topic of this paper. It will be shown that the dependence of 
the S regulator, which controls the rapidity divergence, cancels in the collinear sector. The 
mixing between the UV and IR divergences disappears when the virtual and real corrections 
are added in the collinear sector and in the soft sector. It also disappears when the collinear 
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and soft parts are added in the virtual corrections and in the real gluon emissions. This 
observation implies the separation of divergences to all orders. 

The TMDPDF can be expressed in terms of the integrated PDF (f) q m(z) as 

f q/N (x,k 2 ± ,n) = I — J m (z,\£±,fi)<f> q /N(-,tt) = I —Jm(-,k±,fi)<l> q / N (z,ii). (10) 

■ ' x Z K Z ' Jx Z y Z 



This corresponds to the matching of the operator for the TMDPDF at the scale ~ QX to 
the operator for the integrated PDF at the scale ~ QX 2 , and the kernel J qg (z, k^,//) is the 
matching coefficient. The differential cross section in Eq. ([7]) can be expressed in terms of 
the integrated PDF as 



da 
d 2 q± 



= J dx 1 dx 2 a H(Q 2 , /j,) J d 2 k 1± d 2 k 2± d 2 'q ± 5^ 2 \k 1± + k 2± - rj±- q±)S(rf± ) //) (11) 

Jx\ Z\ Jx 2 Z 2 x Z\ ' v Z 2 ' 



= f dxidx 2 a H(Q 2 ,fi) f — f —W(z 1 ,z 2 ,q ± ,fi)<p q/Nl ( — ,fi)(pg /N J — ,fi). 
J J xi Zi Jx 2 z 2 ^Zt 1 K z 2 > 

where the TMD kernel W(zi, z 2 , q±, ji) is defined as 



W(z 1 ,z 2 ,q±,ii) = J d 2 k 1± d 2 k 2± d 2 r) ± 5 ( - 2) (k 1± + k 2± - rf±- q±) 

xS'(r7_ b n)J qq (zi, ki ± , /j,)J qq (z 2 , k 2± , //). (12) 

All the information on the transverse momentum resides in W(zi, z 2 , q_i_), which should be 
IR finite. It will be computed at next-to-leading order. 
The anomalous dimension for W at one loop is given as 

d 

lw = /J,- r W(z 1 ,z 2 , q±, /i) 

= 7s (qi)*(l - zi)S(l ~ z 2 ) + lj qq {z x ,ql)8{l - *a) + 7j„(*2, ql)*(l - (13) 

where 

7s = /i^5'(q J .,/x) > lj qq = ^J qq {zi,q±,li), U m = n ^Jqq{z 2 , q±, y). (14) 

The anomalous dimensions 75-, 7j 99 and jj contain IR divergences, hence not physically 
meaningful. But the sum is independent of the IR divergence. 

Summarizing our approach, all the particles are on their lightcones, and we define a 
universal TMDPDF in terms of collinear fields only, and employ the soft zero-bin subtraction 
in the TMDPDF to avoid double counting. It is different from other approaches in which 
soft parts are involved in the definition of the TMDPDF. If soft parts are included in 
the TMDPDF, its definition cannot be universal since soft parts are different in different 
scattering processes in general. 
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III. RADIATIVE CORRECTIONS OF THE TMDPDF 



In SCET, the TMDPDF at the parton level is defined as 



fq/N( ,kjj 



P 



(N(p-)\x n ^5(u - n ■ V)5^(k ± - V ± ) Xn \N(p-)). 



(15) 



It is implied that u> > 0, and the longitudinal momentum fraction x is given by x = u>/p~ . 
The TMDPDF fq/N can be treated in a similar way. As emphasized, it consists of the 
collinear fields only, which enables a universal definition of the TMDPDF. There is rapidity 
divergence, and it is regulated by the 5 regulator in the collinear Wilson lines. After the 
zero-bin subtraction, the radiative correction is free of rapidity divergence as in the full 
theory. 

The integrated PDF can be obtained by integrating the TMDPDF over the transverse 
momentum to all orders in a s as 



X 



V- 



(N(p-)\x n ^S(u - n ■ V) Xn \N(p-)). (16) 



This is obviously true at tree level since the tree-level PDFs are given by 

f%(x,kl) = S(l-x)8^(k ± ), <$ N (x) = 8(l-x). (17) 
But one of the main issues involved in TMDPDF is whether it holds true at higher orders 



We are going to construct a formalism in which Eq. (16) works to all orders, and show 
the explicit result at one loop. 



The 8 regulator is introduced in the collinear Wilson lines W n and Wn as 
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perm 



n ■ V + h 



n ■ A r 



ex p 
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perm 



n ■ V + 5 2 



-n ■ A~ 



(18) 



The 8 regulators, as discussed, are employed to regulate the rapidity divergence in TMD- 
PDF. These regulators can be either arbitrary parameters independent of kinematics, or 
the quantities depending on the details of the jets other than the n or n directions. In the 
special case for a back-to-back current, if we put n and n-collinear particles slightly off shell 
by p\ and p|, they are given by 5i = p\jn ■ p 2 and 5 2 = p\/n ■ pi in the latter case. 

The Feynman diagrams for the radiative corrections at one loop are shown in Fig. [TJ 
Fig. [I] (a) gives 



M a 



5(l-x)5(k] 



-(1 + hA+ln^ 2 



V 



-p z 



In' 



P~ 



+ ln 



In 



Si 



+ 2 



7T" 



-p' p 



■ (19) 



Here we use the relation 5^ 2 \k±) = 5(kj_)/7r since the remaining function is independent of 
the azimuthal angle. Fig. [I] (b) and (c) yield 



M h 



.r 



2tt 2 



7T _ 



+ - In 

6 2 p 



Sy 



X) + 



6(1 -x 



8(1 - x) In — 
P 

dn(l — x\ 



In 



—p A 



x\ax 



l — x)'+ (1 - x). 
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r 1 " 




+ 











M c 



a s C F 
2tt 2 



x 



1-x) 



6(1 



x) In—)}, 



:i-^){5(ki)[ln^-l 



lnx(l 





r 1 I 


+ 







(20) 



In Eq. (20), the /1 2 distribution function is introduced. For a well-behaved test function 



/(k 2 ), the integration involving the /x 2 distribution is defined as 



dkl[g(kl)l 2 f(k] 



dk\g{kl)f{kl 



dk 2 ± g(kl)f(0), 



(21) 



where At is a finite upper limit for the relevant physical processes in consideration. Here 
g(k\) is a function which diverges at k^ = 0. At one loop, there are two types of g(k\) 
given by 

1 , ln(kj//i 2 ) 

: ^2"' 92{K-±) - To • 



(22) 



Of course, any definition of the distribution function will do with an arbitrary upper limit 
in the second integral in Eq. (21) instead of /1 2 as long as it removes the IR singularity at 
k^ = 0. The choice of the /1 2 distribution is somewhat cosmetic in the sense that the n 
dependence of the finite part is the same as the radiative corrections of the integrated PDF. 
And there appears no explicit dependence on Ay. 

Since we require that kj_ ~ 0(QX), the transverse momentum should remain finite. 
It is convenient to obtain the information about the dependence of the TMDPDF on the 
small transverse momentum from experiments. But in the effective theory, it can reach 
infinity and the UV divergence can be extracted. However, it becomes subtle if we allow 
the transverse momentum to have arbitrary values extending to infinity. In this case the 
distribution should have a different form and the UV divergence should be included. This is 
relevant when we try to obtain the integrated PDF from the TMDPDF by integrating the 
transverse momentum from to infinity. Then the definition of the distribution function 
includes the UV divergence, and it is referred to as the "infinity distribution" . This will be 
discussed in the next section when we compare the radiative corrections of the TMDPDF 





(a) 





FIG. 1. Feynman diagrams for one- loop corrections of the TMDPDF (a) virtual corrections and 
(b), (c) real gluon emission. The mirror images of (a) and (b) should be included. 



10 



and the integrated PDF. The properties of these distribution functions are described in 
detail in Appendix A. 

In computing collinear matrix elements, the kinematic region with soft momentum is also 
included. Since the soft contribution is computed separately, the contribution from the soft 
region to the collinear part should be subtracted to avoid double counting. This is referred 
to as the zero- bin subtraction [27] . However, there are two possible kinematic regions for 
the zero-bin contribution: soft particles with momentum p s ~ QX, and usoft particles with 
momentum p us ~ QX 2 . The need for the distinction between the soft and usoft zero-bin 
contributions arises depending on whether the processes in consideration are sensitive to 
soft momentum or not. When the zero-bin subtraction was first applied to a heavy-to-light 
current J27J, the soft interaction is decoupled from the current because the interaction puts 
a heavy or a collinear particle off the mass shell. The resultant effective theory is SCETi. 
In order to go down to SCETn, where the fields have fluctuations of order QX 2 , all the soft 
scales of order QX are rescaled to usoft scales of order QX 2 . Therefore there is no need to 
distinguish soft and usoft zero-bin subtractions because there is a smooth transition from 
soft to usoft contributions as we scale down from QX to QX 2 . 

When there is no smooth transition from soft to usoft region as in the Drell-Yan process 
with transverse momentum of order QX, care should be taken about which zero-bin subtrac- 
tion should be performed. If we scale down from soft to usoft scales, while the transverse 
momentum is fixed qx ~ QX, the resultant usoft interactions cannot produce transverse 
momentum of order QX. In other words, the soft part treats physics with momenta of order 
QX, and a naive computation of the collinear part includes the soft part which must belong 
to the soft sector. Therefore the soft zero-bin subtraction is appropriate to correctly avoid 
double counting in the collinear part. As we mentioned earlier, the usoft Wilson lines cancel 
and there is no usoft contribution. 

The zero-bin contribution for M a is given by 

M a = -7T^-°( l ~ x )°{*-±) o ln TT ~ n ln 97 T ■ ( 23 ) 

The same result is obtained whether the soft or usoft zero-bin contribution is considered 
since there is no distinction in the virtual correction. The major difference comes from real 
gluon emission. The soft zero-bin contribution for M&, responsible for real gluon emission, 
is given as 

^(l-,){5(k ± )(-ln —^ + ^ 

,2 X . . 1 . _l„ U /,,2 



A C» = ^(i--)M)(im 2 ^ + 

On the other hand, the corresponding usoft zero-bin contribution is given by 

^.usoft = - x)5(k ± ) - + - In — — + - In — — + — . (25) 



2tt 2 



e 2 e — p 2 5i 2 —p 2 8\ 



This is exactly the same as M° with the opposite sign. The zero-bin contribution for M c 



is suppressed by A compared to M c and is neglected here. Combining Eqs. (19), (20), (23) 
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and (24), the collinear contributions after the soft zero-bin subtraction are given as 

,2 
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(26) 



Just for reference, M& in the usoft zero-bin subtraction is given as 



M, 



b, usoft 



M b - M, 
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b, usoft 
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~2k 2 



{5(ki)[5(l-x)(--ln 2 
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—p L 
ln(l 
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' 12 

x) 



l-z 
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r 1 
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p- J 



(27) 



Each Feynman diagram for the collinear contribution depends on the 5 regulator, but 
after the soft zero-bin subtraction, the collinear contribution is independent of 8\. On the 
other hand, after the usoft zero-bin subtraction, there still remains the dependence on Si. 
This reflects the fact that there is a mismatch in counting the degrees of freedom in the usoft 
zero-bin contribution and in the soft contribution. Furthermore, note that M° + M° usoft = 



from Eqs. (23) and (25), which confirms that there are no usoft contributions in this process. 



To summarize, the soft zero-bin subtraction is appropriate not because it makes the result 
free of Si, but because there is a definite physical reason. Since the collinear contribution 
with the soft zero-bin subtraction is free of the S regulator, there is no need to relate the S 
regulators in the collinear Wilson lines to the offshellness appearing in the soft Wilson lines. 
This is important because the treatment of the S regulator can be applied to processes other 
than Drell-Yan process, in which there are collinear particles other than n and n directions. 
The collinear Wilson line is obtained by integrating out the degrees of freedom of order Q 
when collinear gluons are emitted from other particles. For example, if there are distinct jets 
or heavy colored particles in the final state, the S regulator for a specific collinear Wilson 
line depends, in general, on the offshellness of other jets. The leading eikonalized form is 
the same irrespective of the source of the emitted collinear gluons. If we try to incorporate 
the S regulator by invoking the offshellness of other particles, the S regulator depends in a 
complicated way on the offshellness of other particles. In our scheme, it is not necessary for 
the S regulator in the collinear Wilson line to have a definite relationship with the offshellness 
in the soft part. The dependence on the S regulator in the collinear part disappears when 
appropriate zero-bin contributions are subtracted, not from the soft part. 

In fact, we do not have to relate the S regulators to the offshellness of the particles. The S 
regulators may be used also in the soft Wilson lines, but we stress that the S regulators in the 
collinear and soft Wilson lines are independent regulators. They can be related only when 
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we obtain the 5 regulators from the offshellness in the back-to-back current. Otherwise, the 
attempt to assign some relations between them may simplify the computation, but it evades 
the consistent treatment of the rapidity divergence. 

The total collinear contribution including the self-energy corrections with the soft zero-bin 
subtraction is given at one loop as 

5(ki) 



M col = 2(M a - M° a ) + 2(M b - M°) + M C + (Z ( + R ( )5(l - x)- 



TV 



a s C F 
2tt 2 
1 + x 2 



+ 



(ln^ 5 -lnx)-(l + ^)^ 



(l-x). 



ln(l — x) 



2 2 



/ 1 + x 2 „ . , , u 
+ ( _^ - 2(5(1 -x) In ' 



— p z 



2(S(1 _,)[_!„ _Ay, ( 28 ) 



where and i?^ are the wave function renormalization and the residue for the collinear 
fermion £ at order a s . They are given by 



Z e = 



a s Cp 1 
4ir e 



, i?£ — 



47T 



(29) 



The collinear contribution, or the TMDPDF is independent of 5. That is, it is free of rapidity 
divergence. It is also true in the radiative corrections for the TMDPDF fq/ N (x,k 2 ± ). Note 
that the mixing between the UV and IR divergences still remains. The mixing is cancelled 
if the soft part is added. 



IV. SOFT FUNCTION AND ITS ONE-LOOP CORRECTIONS 



The soft function is defined as 

SM = ^(0\tT[Y r lY w 5^(rj ± +tV ± )Y^Y n }\0}. 

The soft Wilson lines with the offshellness to regulate IR divergence are given by 



Y n= ex P 



1 



perm 



n-V + Ai +i0 



gn ■ A, 



Y- 

5 1 n 



ex p 



i 



perm 



n ■ V + A 2 + iO 



gn ■ A, 



(30) 



(31) 



The operators n ■ V and n ■ V extract the n- and n-components of the momenta from the 
soft gluon. The appearance of Ai and A 2 in Y n and Yn resembles the 5 regulators in the 
collinear Wilson lines, but they are from different origins. Recall that the soft Wilson lines 
are obtained by integrating out the offshellness of the intermediate states when soft gluons 
are emitted from a collinear particle. Therefore, unlike the collinear Wilson lines, Aj's 
depend solely on the offshellness of the relevant collinear particle, not other jets. And it 
should be noted that these regulators in the soft Wilson lines regulate the IR divergence, 
not the rapidity divergence. We can write A x = pf/pi, A 2 = v\lvt- Only in the case of the 



13 




FIG. 2. Feynman diagrams for soft contributions (a) virtual corrections and (b) real gluon emission. 

back-to-back current, the 8 regulators in the collinear Wilson line are given by A x = 5 2 and 
A 2 = S x . 

The Feynman diagrams for the radiative correction of S(t]±) at one loop are shown in 
Fig. [2] The virtual correction in Fig. [2] (a), and the real gluon emission in Fig. [2] (b) with 
their hermitian conjugates are given by 
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and the total soft contribution is given by 



(32) 



M soft = Ml + Ml 



e AiA 2 
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(33) 



The soft function depends on the IR regulators Aj. Especially, there is mixing between the 



UV and IR divergences. When it is combined with the collinear contributions in Eq. (28) 
with p 2 = p 2 , and another collinear contribution from the TMDPDF in the n direction with 
P 2 — V21 the mixing of the UV and IR divergences cancels out. 



V. TMDPDF AND INTEGRATED PDF 

The integrated PDF is defined also in terms of the collinear fields only as 

hM—) = (N(p~)\Xn^(u-n-V) X n\N(p-)). (34) 

P- z 

Compared to the definition of the TMDPDF, the only difference is the absence of the two- 
dimensional delta function on the transverse momentum. That is, the integrated PDF can 
be obtained from TMDPDF by integrating over the transverse momentum. It looks trivial 
to obtain the relation between the TMDPDF and the integrated PDF, but it is not trivial 
to verify the relation at higher orders [Hj. The complication arises due to the regularization 
methods employed to tame the UV, IR and rapidity divergences. The UV divergence can be 
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typically treated in dimensional regularization. The rapidity divergence is handled by the S 
regulator in the collinear Wilson line with the soft zero-bin subtraction. The IR divergence 
is taken care of only when the collinear and soft parts are summed. In this procedure, the 
collinear and soft parts are entangled because the regularization methods complicate the 
transparent relation between the TMDPDF and the integrated PDF, and it needs special 
care to see the relation at higher orders. 

The important issue in establishing the relation is the range of the transverse momentum. 
When the transverse momentum is to be integrated, the region of integration should be 
determined beforehand. If we confine ourselves to finite transverse momentum, there is no 
UV divergence in real gluon emissions. Since the Drell-Yan process with the transverse 
momentum qx ~ QX is considered, we can limit the region of validity accordingly. In that 
case, we use the idea of /i 2 distribution which is suggested exactly for this purpose. On 
the other hand, the transverse momentum can extend to infinity in the effective theory, 
and the /i 2 distribution is not valid any more. There are two possible ways to establish 
the relation between the TMDPDF and the integrated PDF. First, the radiative corrections 
for the integrated PDF are performed in the usual way, and we modify the /i 2 distribution 
in the TMDPDF to include the infinite transverse momentum. Secondly, in the spirit of 
computing TMDPDF assuming that the transverse momentum never reaches infinity but 
some finite scale A?, the radiative corrections of the integrated PDF are modified such that 
the integration over the transverse momentum is limited to At in real gluon emissions, while 
the virtual corrections are performed as usual. 

Let us first consider the case in which the transverse momentum can go to infinity. The 
radiative corrections to the integrated PDF can be computed using the Feynman diagrams 
in Fig. [T] without the two-dimensional delta function for the transverse momentum. We 
organize the collinear matrix elements for integrated PDF including the zero-bin subtractions 
in terms of the virtual and real gluon emissions as 
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And the total collinear contributions at one loop are given as 
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The result is independent of the 5 regulators, which means that there is no rapidity diver- 
gence in the integrated PDF. 

For the TMDPDF, we need another distribution called the "infinity" distribution, in 
which the finite scale in defining the fi 2 distribution extends to infinity. The function of 
the form 



<7i(ki,e,A) 



2c 



(37) 



(ki + A) 1 ^' 

appears in the one-loop correction, where A is a small quantity. For a regular test function 
/(k 2 ), the integral of /(kj_) multiplied by gi(k 2 ,e, A) is written as 



P oo /*oo P oo 

j( rfki^(ki, e ,A)/(ki) = y o dki 5l (ki,e,0)(/(ki)-/(0))+/(0)y o rfki^(ki, e ,A) 

9l (kl, e, 0)] /(ki) + /(0) (- + In ^) . (38) 
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■6 A- 

In the infinity distribution function, we put A = since the IR divergence is already 
regulated. Compared to the /z 2 distribution function, the only difference is the appearance of 
the pole 1/e. This comes from the UV region, which cannot be reached in the /i 2 distribution 
function. Therefore, (?i(k 2 L , e, A) can be written as 

t' 

■e A - 

Similarly, there is another function expressed in terms of the infinity distribution functions 

as 

k 2 
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(40) 



Here two points should be noted. First, the terms proportional to 5(k^) in the fi 2 distribution 
have the same [i dependence as those in the infinity distribution except the UV poles. 
And secondly, the infinity distribution is not exactly a distribution function in a rigorous 
sense. It is because the UV pole cannot be peaked near k\ ~ 0, therefore it should be 
understood that the infinity distribution function is meaningful only after the integral over 
k 2 is performed. This awkward situation occurs since we are going to compare the TMDPDF 
and the integrated PDF, and the integrated PDF is obtained by integrating k^ over all values 
including infinity. 

The total radiative correction for the TMDPDF is given as 
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If we integrate over kj_, the infinity distributions vanish and we obtain 
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(42) 



which is exactly the same as the radiative corrections for the integrated PDF in Eq. (36). 



Secondly, we consider the case in which the transverse momentum has the range < 
k 2 < Ay. Then the radiative corrections for the integrated PDF is handled differently. 
The loop momentum I is defined in D dimensions, and there is a two-dimensional delta 
function. We introduce the momentum vector in D — 4 = — 2e dimensions, and lj_ is a 
two-dimensional momentum vector. Then I 2 can be written as I 2 = — lj_ — l\. First 
the two-dimensional transverse momentum is fixed in real gluon emissions and the integral 
over l± in D — 4 dimensions is performed first. Then we integrate over the two-dimensional 
transverse momentum. Since the virtual corrections are not affected, M a is the same. And 
the matrix elements for the real gluon emissions are modified to be 
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The whole collinear part of the integrated PDF at one loop is given by 

7 ^2 q ,,2 
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If we integrate Eq. (28) from to A 2 ^, Eq. (44) is obtained. In either case, the integrated 
PDF is consistently obtained by integrating the TMDPDF over the transverse momentum 
to one loop. 



VI. DIFFERENTIAL SCATTERING CROSS SECTION AT NLO 



The TMDPDF can be expressed in terms of the convolution of the integrated PDF and 
the TMD kernel which contains all the information on the transverse momentum in the 



TMDPDF. Its relation, Eq. (10), is given again as 



f dz /X \ f dz /X \ 

f q / N (x,k±,fj) = —J qq (z,k±,fj)<f> q / N (-,ij,)= —J q J-,kx,n)<f> q / N (z,n) } (45) 

Jx Z K Z ' Jx Z y Z ' 

where <p q /N{z) is the integrated PDF. We will use this relation to obtain the renormalization 



group behavior of the scattering cross section to one loop. Eq. (45) can be regarded as 
matching between two SCET's, where the TMDPDF is the quantity in which p ~ QX and 
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the integrated PDF is the one where p ~ Q\ 2 . And the kernel J qq (z, kj_, //) is the Wilson 
coefficient in matching. At zeroth order in a s , Jf® and (j) q °) N are given by 



4?(^k ± ) = t S {l - z)5(kl), ^/ N (z) = 5(1 - z). 



71 



(46) 



The TMDPDF and the soft functions at order a s are obtained from Eqs. (28) and (33) by 



removing the UV divergent terms. They seem to have the mixed divergent terms, but they 
disappear when all the contributions are added, therefore it is of no concern here. They are 



given as 
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(47) 



Expanding Eq. (45) to first order in a s , we obtain 
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The integrated PDF at order a s is given by 
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and the anomalous dimension for the PDF is given by 
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Using Eqs. (48) and (49), the new result is the TMD kernel at order a s , which is given as 
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The factorized differential scattering cross section is written as 
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The anomalous dimension of at one loop is given as 

d 



-]}■ 



lw(zi,z 2 ,q ± ,/j,) 



d 



din fx 

a s C F 5(q 2 ± ) 

2lT 7T 



W(z 1 ,z 2 ,q ± ,Li) = 2li —W(z 1 ,z 2 ,q ± ,n) 



dfi 2 



(54) 



(55) 



2\n^5(l- Zl )5(l-z 2 ) 



Q ^ (qi) [G F (2h 1 ^ +3)<5(l-c 1 )*(l-^) 



2"7T 7T 



P OT (^l)<J(l " 22) " Pqqi^W ~ Zl ) 

also 

lw(zi,Z2,q^_,fi) = 5(q 2 ± ) ; y w (z 1 ,z 2 ,Li) 



where P qq (z) = Pqq(z). It can be also written as 



(56) 



= l.J qq (zi,q 2 ± ,fi)5(l - z x ) +1j„(z 2 ,<£ l ,h)5(\ - z 2 ) 
+ls (ql)5(l- Zl )5(l-z 2 ). 

Here ^j qq1 7j 5(? and 75 are the anomalous dimensions for J qq , J qq and the soft function, which 



Jqq J I -Jqq 

are given as 



-yj(x,q 2 ± ,n) 



d j / N a s C F 5(q 2 ± ) 

J(x,q ± ,n)~ 



dlnfj, 



2tt 



^[25(1 -x) In 



—p 2 



1 + x 2 



(1-x). 



, 2 \ d Q ( \ n a sC F 5(q 2 L ) Li 2 



din// 



7T 



7T 



AiA 2 ' 



(57) 



19 



where Ai = p\/n ■ p±, and A 2 = p 2 /n ' P2, and J = J qq = Jqq. Note that jj and 75 still 
depend on the IR-dependent terms, but the sum jw has no IR divergence. Compared to 
the anomalous dimensions corresponding to 7j and 75 in Ref. |16j . the dependence on the 
rapidity scale v is tantamount to the IR dependence in our approach. 
The hard coefficient C(Q,fi) is given by 

C(Q, ,) = 1 + ^ (- in 2 ^ - 3 In - J - 8 + £) , (58) 

and the anomalous dimension of the hard part H(Q,p) = \C(Q,fi)\ 2 is given as 

7 „M = -^(21„! + 3). (59) 



Combining 714/, 7# and 7^, it is clear that the factorized scattering cross section, Eq. (53), 
is independent of the renormalization scale fi at one loop. The numerical analysis for the 
scattering cross sections is interesting with the solution for the renormalization group equa- 
tion. This can be applied to the processes with small transverse momentum such as Drell-Yan 
process, Z boson production as well as Higgs production both at Tevatron and LHC. The 
numerical estimates involving TMDPDF in these processes will be presented in future pub- 
lication. 



VII. DECOUPLING OF THE UV AND IR DIVERGENCES 

In order to see how the decoupling of the UV and IR divergences is achieved in various 
combinations, it is convenient to collect all the results obtained so far. The results for the 
integrated PDF are summarized as follows: The collinear parts are separated into the virtual 
and real parts as 



'PDF l 



5(1 



I R 

'PDF 



(Pt) 



2tt 

u s Cf 
2tt 
1 + x 2 



X 



2 3 2, u 2 
+ — + - In ■ 



e 2 ' 2e ' e "* —p\ 



, 3 i ^ 



7 



2 u 2 



S(l-x){- 



-In 



2 ^ 



In x + In — 



e -p\ 
2 



-Px 



+ (l + x 2 



7T Z 

~ ~6 
ln(l 



-Pi 
2(1- 



Ti 

~6 



x. 



1 



X 



(60) 



Note that there is another collinear contribution from the n direction, which is obtained 
from replacing p\ by p\. Similarly, the soft function is given as 
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The mixing terms (ln/z 2 /(— p 2 ))/e exists in Mp DF , Mp DF , M^ ft and M^ ft respectively. How- 
ever, if we add the collinear parts Mp DF + Mp DF , the mixing term cancels. It is also true for 
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the soft parts M s ^ ft + M s ^ ft . Therefore the total contribution is free of the mixing. There is 
another combination in which the mixing cancels. They are Mp DF (p1) + Mp DF (pl) + M s ^ ft , 
and Mp BF (pf) + Mp DF (pl) + M^ ft . It means that the virtual corrections of the collinear 
and soft parts do not contain the mixing of the UV and IR divergences, nor the real gluon 
emissions of the collinear and soft parts. 

In the case of the TMDPDF, the virtual corrections of the collinear and soft parts have 
the same expression except 5(kj_), hence the mixing is not there. In the real gluon emissions, 
the possible mixing terms appear when approaches infinity (or A T ) in the terms 



u 2 u 2 1 

M tmdpdf = 2 Mb + M C ^ -2 (in + In 



-pi -P2 

,2 ,,2 



k 



< ft = ^2(ln^ + ln^)[^ 2) (62) 



where Ai = p\/n ■ p\ and A2 = p\jn ■ p2- Other terms may contain divergences, but the 
UV and IR divergences are decoupled. As can be seen clearly, these are cancelled when the 
real gluon emissions in the collinear and soft parts are added. The mixing of the UV and 
IR divergences disappears in the virtual corrections and in the real gluon emissions of the 
collinear and soft parts respectively. 

This observation is useful since the virtual corrections of the TMDPDF have the same 
form as those of the integrated PDF except 5(kj_). And they do not have the mixing term. 
Therefore the real gluon emissions should not have any mixing term either, when integrated 



over the transverse momentum. As can be seen in Eq. (62), the mixing term is cancelled 
before the integration over the transverse momentum. If it is true to all orders in a s , we can 
safely use this fact in order to obtain the anomalous dimension of the combined collinear 
and soft parts. 



VIII. CONCLUSIONS 

The divergence structure in the TMDPDF is intricate because the rapidity divergence 
should be also handled as well as the usual UV and IR divergences. One of the major 
problems in studying TMDPDF is to find a regularization method to treat the rapidity 
divergence. The dimensional regularization can handle the UV divergence, and possibly 
the IR divergence, but it cannot isolate the rapidity divergence. As far as the UV and 
IR divergences are concerned, it is convenient to use dimensional regularization to extract 
the UV divergence, and to employ the offshellness of external particles. There are various 
ways to treat rapidity divergence, and we have found that a simple way to regulate the 
rapidity divergence is achieved by inserting the 5 regulators in the collinear Wilson lines only. 
Individual collinear Feynman diagrams depend on the regulator 5, but the total collinear 
contribution with the zero-bin subtraction is independent of the regulator, hence no rapidity 
divergence. The regulators similar to 5 may be inserted, say, in the soft Wilson lines, or 
in the collinear fermion propagator, but they regulate the IR divergence, not the rapidity 
divergence. 
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For comparison, let us recapitulate previous approaches in handling the rapidity diver- 
gence. Collins and Soper [T2] tilted the lightcone such that the rapidity divergence is handled 
by an extra scale in association with this tilting. Chiu et al. [16J modified the collinear and 
soft Wilson lines by introducing a rapidity scale u, which treats the rapidity divergence. In 
both approaches, the hard, collinear and soft parts depend on the additional scale, while 
the cross sections do not. Therefore another renormalization group equation with respect to 
this scale is introduced. Another approach is to choose the axial gauge [281 I2H] , for example 
n ■ A n = 0, instead of covariant gauges such as the Feynman gauge we chose. The advantage 
of this gauge choice is that many Feynman diagrams from the Wilson lines vanish, but there 
should be additional transverse gauge link which links the operators for TMDPDF at infinity 
so that it becomes gauge invariant. Ref. [22] treats the rapidity divergence in a similar way 
to ours, but they insert the 5 regulators in both the collinear and the soft Wilson lines and 
give a specific relation between the collinear and the soft regulators. That relation, as we 
point out, holds only for the back-to-back current. 

The separation of the UV and IR divergences is achieved in the sum of the virtual and real 
corrections in the collinear sector and in the soft sector, and in the sum of the collinear and 
soft parts in the virtual corrections and in the real gluon emissions. In effective theories like 
SCET, the mixing shows up often in intermediate steps, but it cancels in the final result. 
In the radiative correction of the TMDPDF, it also holds in the sums mentioned above. 
Therefore the renormalization group scaling can be considered only for the hard part and 
the sum of the collinear and soft parts. This prescription seems to violate the factorization in 
its strictest sense that the collinear and the soft parts are separately well defined. We rather 
claim that the attempt to combine the collinear part with part of the soft part such that 
the UV and IR divergences are decoupled is futile in its strictest sense that no such thing 
occurs with general regulators unless a very specific set of relations is imposed. The soft part 
interacts with all the collinear sectors in the process, hence contains the information of all 
the collinear sectors through their offshellness. Therefore it is physically more reasonable to 
consider the scaling behavior of the sum of the collinear and soft parts in order to decouple 
the UV and IR divergences. This may sound cumbersome, since it means that we have 
to compute all the collinear and soft parts first in complicated processes with multijets or 
heavy colored particles. But we stress that it is the correct procedure to treat high-energy 
scattering. It will be interesting to see if our procedure also works in other processes involving 
small transverse momentum. Various aspects of our formalism will be investigated in future 
publications. 

The important result in this work is the factorization formula for the differential scattering 
cross section in Drell-Yan process with small transverse momentum, given by Eq. (11 ). The 
scattering cross section is factorized into the hard part, the TMD kernel and the PDF's. 
The hard part is the matching coefficient between the full theory and the SCET with the 
scale QX. The TMD kernel is obtained by matching the two SCET's with the scales QX 
and QX 2 . Therefore the TMD kernel should be free from IR divergence to all orders, and 
it is shown explicitly at NLO. Note that the TMD kernel W(z\, Z2, q_i_, fj) is expressed in 
momentum space, rather than in impact parameter space. At least at NLO, the anomalous 
dimension of W is proportional to 5(qi), hence the dependence of the scaling behavior on 
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the transverse momentum is trivial. 
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Appendix: ji 2 distribution functions and infinity distribution functions 



Suppose that a function g(x) diverges as x — > 0. Then the plus distribution function 
can be defined, in general, as 

a pa rb 

[g(x)] b f(x)dx = / dxg(x)f(x) - /(0) / dxg(x) 
o Jo Jo 

9(x)\f(x) - /(0)j - /(0) f b dxg(x), (A.l) 

for a regular test function f(x). The need for the plus distribution function is to separate 
the IR divergent part. When g(x) diverges as x — > 0, the integral is tamed to be finite 
by subtracting /(0). In this general definition of the plus distribution function, the upper 



limits a and b are arbitrary as long as the lower limits of the two integrals in Eq. (A.l) are 
both zero to guarantee that the IR divergence at x = is cancelled in this subtraction. The 



finite remainder, the second term in Eq. (A.l), depends on the limits a and b. In a special 
case a = b = 1, the familiar plus distribution function is obtained as 



dxg{x) f(x) - f(0) 



(A.2) 



There are two types of functions in the radiative corrections, which are of the form 

k 2 



gi(k 2 ± ,e,S) 



2c 



(ki + A)^ 



g 2 (k 2 ± ,e,5) 



/i^ln 



A 



(ki - A)i- 



(A.3) 



where A approaches zero. If A = 0, these functions are singular at ki = 0, that is, IR 
divergent. We extract the IR divergence as the coefficient of 5(ki), and the remainder is 
the distribution function free of the divergence. With the presence of the nonzero A, the 
IR divergence is regulated by A instead. However, there is a big difference in considering 
the distribution functions g\ and g-i in ki-space compared to [g(x)]6 since the upper limit 
can reach infinity. Irrespective of whether the upper limit goes to infinity or not, g\ and 
#2 diverge at ki = with A = 0. If the upper limit in considering the integral including 
gi or <?2 reaches infinity, the integral also has an UV divergence. Therefore we classify 
the distribution functions according to the upper limits. These two cases have their own 
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physical importance. If we restrict the transverse momentum to be of order QX and collect 
experimental data, there is finite upper limit which we call At- The relevant distribution 
function is called the "/x 2 distribution". In the effective theory, can reach infinity and 
the UV divergence from the integration offers the information on the scaling behavior. The 
"infinity distribution" is devised for this case. 

First let us consider the case in which the transverse momentum remains finite up to 
some scale A T - If a regular test function /(k 2 ) multiplied by gi(k 2 , e, A) is integrated from 
to A T , it yields 



/ T rfki^(ki, e ,A)/(ki) 
Jo 

2 2 

= / T rfki^(ki,6,A)/(ki)- f dkfoi(ki,c,A)/(0)+ f dki 5l (k 2 ± ,e,A)/(0) 

./ii ./ii Jo 
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9l (k 2 ± ,0,0) J(ki) + /(0) In 



(A.4) 



In the final expression, we put e = 0, and A = in g(k 2 ± , e, A) because the integral has no 
pole in e, and the possible infrared divergence near k 2 ~ is cancelled in the subtraction. 
This defines the fj, 2 distribution. g 2 can be treated in a similar way to yield 



<?i(k 2 ,e,A) = 
(? 2 (k 2 ,e,A) = 
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(A.5) 



Note that these functions are independent of \x since we can put e = from the beginning. 
However, the \i dependence is split into the delta function and the ji 2 distribution function 
for convenience. The independence can be seen by computing the integral 
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(A.6) 



If the upper limit for the transverse momentum reaches infinity, g± and g 2 contain UV 
divergences: 

,2 



roc I i, z 

/ dklg l (kl,e,A) = - + \n fJ - + 0(e) 
Jo e A 
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Therefore g\ and g-i can be expressed as 



<?i(ki,e,A) 
</ 2 (kl,e,A) 



In 



A 



5(kl) 
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<7i(kl,e,0) 



(A.8) 
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Note that Eq. ( A.8 ) is not strictly correct since the UV divergence is not peaked near k]_ ~ 0. 
Therefore they are meaningful only inside the integral over the transverse momentum which 
reaches infinity. Then the UV divergence is correctly accounted for. 
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